Backward waves in a grounded bilayer slab containing 
double-negative (DNG) and double-positive (DPS) metamaterials 

E. Cojocaru 
Department of Theoretical Physics 
Horia Hulubei National Institute of Physics and Nuclear Engineering 
P.O.B. MG-6, 077125 Magurele, Romam^ 

Abstract 

Simple dispersion relations for the guided modes in a grounded DNG/DPS bilayer slab are given 
in terms of normalized parameters. Relations corresponding to the grounded single-layer DNG 
slab are refound as specific cases. Numerical examples are given showing dispersion curves of the 
lower order modes and the respective total normalized power carried on the propagation direction. 
Snapshots obtained by the finite-difference time-domain method are provided showing the electro- 
magnetic field inside the grounded DNG/DPS bilayer slabs. Since an important characteristic of 
the guided modes in the slab containing a DNG layer is the existence of a turning point (TP) at 
which the power carried by each mode of order m > equals zero and changes the sign, we present 
implicit relations at the TP for the normalized parameters of the guided modes in the grounded 
DNG/DPS and DNG slabs. We show that a thin DPS layer coating on the grounded DNG slab 
produces a shift of the TP on the dispersion curve. 
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I. INTRODUCTION 



Recently, the materials having negative permittivity and negative permeability, which 
are labeled as double-negative (DNG) materials, have attracted much attention, especially 
because of their unusual physical properties, which are different from those of conventional 
double-positive (DPS) materials. Veselago jl| was the first to study theoretically the DNG 



materials, but these materials became a topic o 
demonstration of their negative refractive index 



' high interest after the first experimental 
2J. Much attention has been paid to their 



potential applications, a pioneering work being Pendry's suggestion of perfect lenses made 
of DNG materials Characteristics of these materials can be found for example in jj- 
Electromagnetic wave propagation in these materials has been extensively investigated 



theoretically and experimentally 



are capable of supporting 
described for example in 



Theoretical studies have shown that such materials 
3ackward waves jsjj]. Waveguides containing DNG layers are 



12-18 



materials have been analyzed in 



Guided waves in grounded single-layer slabs of DNG 



12,ll9r|22|. 



In this paper we consider grounded slabs containing DNG/DPS bilayers. Characteristics 
of the single-layer DNG slabs are refound as specific cases from those for the bilayer slabs. 
Dispersion relations are given in terms of normalized parameters. Numerical examples show 
dispersion curves of the lower order modes and the total normalized power carried on the 
propagation direction by the respective modes. Examples of electromagnetic fields inside 
the bilayer slab are also given. Snapshots obtained by the finite-difference time-domain 
(FDTD) method 23| are provided. An important characteristic of the guided modes in the 
slab containing a DNG layer is the existence of a turning point (TP) at which the power 
carried by each mode of order m > equals zero and changes the sign. We present implicit 
relations for the normalized parameters of the guided modes in the grounded DNG/DPS 
and DNG slabs at the TP. 



II. GENERAL RELATIONS 

The geometry of the grounded bilayer slab is shown in Fig. [IJ The slab of thickness d\ + d2 
is infinite along the y and z directions. It is grounded by a perfectly electric conducting 
substrate in the plane x — 0. The relative material constants are (e_j, jUj), with j = 1,2, for 



1 



the layers of the slab, (e c , /x c ) for the cover medium, and (e s ,// s ) for the substrate, so that 
£iA*i > £2/^2 > CcyUc- For simplicity, the analysis in this paper is restricted to the case of 
frequency-independent and lossless material constants in order to obtain relations between 
the geometrical and material parameters. Only in the FDTD numerical simulations we use 
the Drude-Lorentz model for the material constants of the DNG layer 23J. We consider 
the modes with a time harmonic variation exp(iu)t) propagating in the z direction, their 
electromagnetic fields varying as expi(ut — (3z), where (3 is the modal phase constant. For 
transversal magnetic (TM) modes, the magnetic field is along the y direction, whereas it is 
the electric field along the y direction in the case of transversal electric (TE) modes. The 
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FIG. 1: Geometry for the grounded bilayer slab of layer thicknesses dj and relative material 
constants (ej,fj,j), where j = 1,2, with a cover medium of relative material constants (6 c ,/i c ). 

magnetic field H y for the TM modes takes the form 



< x < d\ 
di < x < di + d 2 
x > d\ + 0?2 



H y = cos(k tl x), 
H y = Ai cos[k t2 (x - di)] + Bi sm[k t2 (x - di)], 
H y = A 2 exp{-K c [x - (di + d 2 )}}, 



(1) 



where the factor exp(— i(3z) is skipped, and 



Ax = cos(k t xdi), B x = sm(k t idi), A 2 = A x cos(k t2 d 2 ) + B x sm(k t2 d 2 ), (2) 

k t 2t\ 

k t j = (klejfij — /3 2 ) 1 / 2 when (3 2 < ejfij, 

k t j = —i(/3 2 — klej/dj) 1 ^ 2 when (3 2 > ejfij, (3) 

with j = 1,2, j3 = 0/k o is the relative phase constant, k = (^(eo/io) 1 ^ 2 is the wavenumber 
in vacuum, and 



k c = (/3 2 - k 2 Q e c n c ) 



1/2 



(4) 



Similarly, one can express the electric field E y for the TE modes. The total normalized 
power carried by each mode on the z direction takes the form 



P = (Pi + P 2 + P c )/(|Pl| + \P 2 \ + \Pc\), 



where Pi, P 2 , and P c are given up to a common factor by 



Pi oc — 
ei Jo 



\H y \ 2 dx, 



1 



d!+d 2 



1 



P2 oc — / \H y \ 2 dx, P c oc — 



e 2 Jdi c c Jd L +d 2 

for the TM modes. Similarly, one can express Pi, P 2 , and P c for the TE modes. 



(5) 



\H y \ 2 dx, (6) 



III. DISPERSION RELATIONS 



Dispersion relations are determined either from the poles of the complex-amplitude re 
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25] taking 



flection coefficients expressed in terms of the transfer matrix elements 
into account that e s = —00 and fi s = 26| . or by applying the boundary conditions of the 
tangential electromagnetic fields at interfaces. One obtains the following dispersion relation 

e c k t i 



tan(k t idi 



for the TM modes, and 



ta.n(k tl di 



ti 



1 - ^^tan( k t2 d 2 ) 
e 2 K c 



1 + ^r^- ta,n(k t2 d 2 ) 



I 



1 + -*r^- ta,n(k t2 d 2 ) 



e c k 



12 



(7) 



12 



I 



lick 



12 



fi 2 K c 



t&n(k t2 d 2 ) - 1 



(8) 



for the TE modes, where k t j, with j = 1,2, and k c are given by relations ([3]) and 
respectively. In the case (3 2 < e 2 fi 2 < ei/ii we introduce the normalized parameters 



u c = (/? 2 - e c /x c ) 



1/2 



u i 



(ejfij - /3 2 ) 1/2 , vj = k dj, j = 1, 2, 



7j = ejuj (e c Uj) for TM modes, r )j = —ii c Uj/(njU c ) for TE modes. (9) 



1 - tdJl{u 2 V 2 ) / ^ 2 



Using notations ([9]) in ((?]) and ([8]) gives 

tan(wit>i) 

71 1 + 72tan(M 2 u 2 ) 

Thus, the TM and TE modes satisfy the normalized dispersion relation 



(10) 



1 

Vi = — 
Ui 



tan" 1 7i[l-tan(^ 2 )/ 72 ] + ^ 
1 + 7 2 tan(-u 2 -u 2 ) 



m = 0,1,2, 



(11) 
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When d 2 = 0, then v 2 = 0, and (TTTT) reduces to the simple form 

v x = (l/«i)(tan- 1 7i + ^), / = 0,1,2,... (12) 

for the TM and TE modes, the cutoffs (3 = (e c /i c ) 1//2 being at v\ = In, with I = OA, 2, . . . , 
for the TM modes, and vx = (21 - l)vr/2, with I = 1, 2, . . . , for the TE modes [lj. Note 
that, for simplicity, we label the modes by I or m without any reference to the standard 
mode labeling of DPS slabs. In the case e 2 /i 2 < f3 2 < e^Ux, tacking into account relation Q, 
we denote 

u 2 = -iu 2 , u 2 = (f3 2 - e 2/ u 2 ) 1/2 , 72 = «72, 
7 2 = e 2 M c /(e c w 2 ) for TM modes, 72 = /A^ 2 /(/i 2 w c ) for TE modes. (13) 



Then, the TM and TE modes satisfy the normalized dispersion relation 



1 

vi = — 
tti 



_ x 7i[l + tanh(M 2 w 2 )/7 2 ] 
tan — — - — r h m7r 

1 + 7 2 tanh(M 2 -u 2 ) 



(14) 



In the case (3 2 > e\fj,i > e 2 /x 2 both 71 and 7 2 have imaginary values, but there are real 
solutions to j3 of the TE mode [18]. We denote u\ and 71 like in ( |T3l) . Then, the evanescent 
mode TE of the grounded DNG/DPS bilayer slab satisfies the normalized dispersion relation 

7i[l + tanh(u2i>2)/72] 



v\ = — tanh 1 

Ui 



(15) 



1 + 7 2 tanh(w 2 f 2 ) 

Since 71 < 0, V\ in f lT5|) has real values, with v\ — > when j3 — > 00. When v 2 = 0, relation 
( JT5l) becomes 

Vl = (l/u 1 )tajih~ 1 (-^ 1 ). (16) 

IV. NUMERICAL EXAMPLES 

For simplicity we consider the cover medium is air (e c = \x c = 1). Dispersion curves are 
shown in Figs.|2^a) and (b) for the TM and TE modes, respectively, when (ei, /ii) = (—4, —2) 
and (e 2 , /z 2 ) = (2, 1.5), at v 2 = 0, 1, and 2. The dispersion curves overlap in the upper region, 



at (3 2 > e 2 /i 2 , but in the lower region, at (3 2 < e 2 \x 2 , there are composed modes [16J. Thus, it 
is more easy to look after the modes starting from the upper region of the dispersion curves, 
at (3 2 > e 2 ji 2 . In Fig. Mjo), the TE mode allows real solutions at f3 2 > ei/ii, the starting 
point being marked by an asterisk. For each mode of m 7^ 0, the total power carried on the 
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FIG. 2: Dispersion curves for (a) TM and (b) TE modes in a grounded DNG/DPS bilayer slab of 
material constants (ei,/zi) = (—4,-2) and (€2,^2) = (2,1.5), when V2 is constant: V2 = (dashed 
line), 1)2 = 1 (thin line, marker .), and v 2 = 2 (thick line, marker +). Only the lower order modes 
are considered: m = (red), m = 1 (blue), m = 2 (magenta), m = 3 (cyan), and m = 4 (green). 
The TP of each mode is marked by a small circle. The total normalized power P carried by 
each mode on the propagation direction is represented in (c) and (d) for the TM and TE modes, 
respectively. 



z direction equals zero at the TP, as shown in Figs. [2](c) and (d). Note that the TPs are 
more distinctly seen on the P against v± curves than on the dispersion curves. We denote 
the values of (3 and Vi at the TP by 0tp and t>i TP , respectively. At (3 < /?tp, the total power 
P is positive (forward wave), whereas at > /3 tp, the total power carried by each mode of 



m 7^ is negative, (backward wave) [ll|, [12|, [17|, [18]. In Fig. [2J one can see that, the greater 
is t> 2 , the greater are /3tp and fi TP . At given v 2 , the greater the order of the mode, the 
smaller 0tp and the greater t>i TP . The fields H y and E y in the slab are shown in Figs. [3]^a) 
and (b) for the TM and TE modes, respectively. At 2 < €2^2 the field is maximum in the 
DPS layer, while at p 1 > €2^2 the field is maximum in the DNG layer. Snapshots of TE 
backward propagating waves obtained by the FDTD method 23[ are shown in Figs. H(a) 



and (b) for the grounded DNG and DNG/DPS slabs, respectively. The Drude-Lorentz 



model was used for the dispersive material constants of the lossless DNG layer 



whereas 



frequency-independent and lossless material constants were considered for the DPS layer. 
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(a) (b) 

FIG. 3: The fields H y in (a) and E y in (b) for the TM and TE modes of a grounded DNG/DPS 
bilayer slab at the points shown in Figs. EJa) and (b) by black markers, when v 2 = 2. For the TM 
mode in (a), v\ = 1.5, with j3 = 1.48 (full line), and j3 = 2.53 (dashed line). For the TE mode in 
(b), vi = 2.5, with p = 1.26 (full line), and $ = 2.4 (dashed line). 

The computation space is 160 x 60 cells and 200 time units is approximately the time to 
run a simulation. A perfectly matched layer was placed at the upper side of the simulation 
domain. The slowly rumped continuous sinusoidal source is placed at the center of the DNG 
layer. We considered vi = 1.8 with u/(2nc) = 0.28 in Fig. Ufa), and v\ = 2.25, d 2 = di/2, 
and u/(27rc) = 0.36 in Fig. IH(b). Contrary to the forward propagating waves in a classical 



sla 

nj. 



d, the backward waves are propagating from the lateral sides toward the centered source 



pi 4tp, ,| 


94 ill? 44 

■iliilti it! 





(a) (b) 

FIG. 4: Snapshots of TE backward propagating waves (a) in the grounded lossless DNG slab of 
frequency-dependent material constants (ei,/ii) = (—4, —2) with v\ = 1.8 and uj/(2ttc) = 0.28, and 
(b) in the grounded lossless DNG/DPS slab of frequency-dependent material constants (ei,/xi) = 
(—4,-2) and frequency-independent material constants (£2,^2) = (2,1.5), with v\ = 2.25, d 2 = 
di/2, and w/(2irc) = 0.36. The sinusoidal continuous source is placed at the center of the DNG 
layer. 



Now, consider a grounded DNG/DPS bilayer slab of relative material constants (ei, /ii) = 
-4, —2) and (e 2 , ^2) = (4, 2). Note that d\ 7^ d 2 - Dispersion curves are shown in Figs. E^a) 
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FIG. 5: Dispersion curves for (a) TM and (b) TE modes in a grounded DNG/DPS bilayer slab of 
material constants (ei, /Ui) = (—4, —2) and (e2, ^2) = (4, 2), when V2 = (dashed line)and V2 = 0.5 
(thick line, marker +). The colors are kept the same like in Fig. [2J The TP of each mode is marked 
by a small circle. The total normalized power P carried by each mode on the propagation direction 
is represented in (c) and (d) for the TM and TE modes, respectively. 

and (b) for the TM and TE modes, respectively, at v 2 = and 0.5. The total normalized 
power P carried by each mode on the propagation direction is represented in (c) and (d) for 
the TM and TE modes, respectively. Dispersion curves at v 2 7^ are translated as against 
those at v 2 = 0, so that, for a given mode of m / 0, 



1)2=0 



V2=0 



+ v 2 . 



(17) 



Now, consider a grounded DNG/DPS bilayer slab of relative material constants (ei, /ii) = 
(—4,-2) and (62,^2) = (2,4). Dispersion curves are shown in Figs. [6]^a) and (b) for the 
TM and TE modes, respectively, at t>2 = 2. The total normalized power P carried by each 
mode on the propagation direction is represented in (c) and (d) for the TM and TE modes, 
respectively. Intricate modes exist at V\ > t>2, as for example, the first mode at v\ > t>2 
has three TPs. The total normalized power P has positive values between the two TPs of 
the same order m, as shown in Figs. [6]^c) and (d). The modes at v\ > t>2 have P < 0, 
whereas the modes at v\ < t>2 have P > 0. In Fig. [6] (d), the mode of m = 2 has P ~ 
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FIG. 6: Dispersion curves for (a) TM and (b) TE modes in a grounded DNG/DPS bilayer slab of 
material constants (ei,y«i) = (—4,-2) and (€2,^2) = (2,4), when V2 = 2. The colors are kept the 
same like in Fig. [2j The TPs are marked by small circles. The total normalized power P carried by 
each mode on the propagation direction is represented in (c) and (d) for the TM and TE modes, 
respectively. 




FIG. 7: Snapshot of the TE backward propagating wave in the grounded lossless DNG/DPS slab 
of frequency-dependent material constants (ei,//i) = (—4, —2) and frequency-independent material 
constants (€2,^2) = (2,4), with v\ = 2, cfe = d\, and lo/(2ttc) = 0.32. The sinusoidal continuous 
source is placed at the center of the DNG layer. 

ue of 



23j is 



at v\ = 4.2282, but P does not change the sign, and so, there is not a TP at that va_ 
v\. Snapshot of the TE backward propagating wave obtained by the FDTD method 
shown in Fig. [TJ when v\= V2 = 2, with cj/(2ttc) = 0.32. After a long time, more than 300 
time units, the field stays concentrated near the continuous centered source. 
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V. IMPLICIT RELATIONS AT THE TP 



Since the TP is an important characteristic of the modes in the grounded DNG/DPS and 
DNG slabs, here we present implicit relations for determining /3tp when i> 2 is constant. 

A. The case (3 2 < €2^2 < eiMl 

At the TP, Vi is minimum and dvi/d(3 = 0. Thus, using ( 1TT|) gives the following implicit 
relation 



7i(! ~ W72) 
1 + 72^2 



tan 



(T71F 



L(l + 72t2)2 + 7 2(l-t2/72) 2 



where jj, with j = 1, 2, is defined by relation (EJ), ^2 = tan(ii 2 i>2), 

a — — 1 for TM modes, a = 1 for TE modes, 



?7ic 



«c " «2 

with u c defined by and 

V12 = ei/ii - e 2 /i 2 , 



, (72 - -) - <7-^(l + t 2 )(72 + _) + -^A^, 



72 



"2 



72 



.7 = 1,2. 



(18) 

(19) 
(20) 

(21) 



Once /?tp is determined with ( fTHj) . the respective value t>i TP is obtained with ( TTT1) . When 
t>2 = 0, the implicit relation (lisp becomes 



71 = tan {a-fir] lc /[u 2 c (l + 7?)]} • 



(22) 



This is the implicit relations for determining /3tp in the case of a grounded single-layer DNG 
slab. 

B. The case e 2 /^2 < /3 2 < eifi\ 



Using dvi/dp = in ( |14p gives the following implicit relation for determining /3tp when 
i>2 is constant, 



7i(! + W72) 
1 + 72*2 



tan 



(l + 7 2 t 2 ) 2 + 7i 2 (l+^/72) i 



where 72 is defined by relation ffTB"]) . t 2 = tanh(M 2 f2), and 



— (1 + t 2 ) + ^o-(72 + — ) + CT-^ 



«2 



72 



1 _ ?)R2 _i_) + 2 "^ 

"2 " 72 



M 2 U2 



Once /9 TP is determined with ( )23|) . the respective value t>i TP is obtained with ([T 



(23) 



(24) 
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C. Numerical examples 



Using implicit relations ( fl8|) and ( 123|) gives at v 2 = 1 in Fig. [2]^a) the following values of 
\Ptp,v 1tp ]: [1.8782,1.1329] at m = 1, [1.5375,2.5163] at m = 2, [1.3741,3.8101] at m = 3, 
and in Fig. 12(b), the following values: [1.5496, 1.9174] at m = 1, [1.3311,3.2039] at m = 3, 
and [1.2266,4.4481] at m = 4. In Fig. [6](a), we obtain for the first intricate mode at 
v x > v 2 the following values of [(3 tp ,v 1tp ]: [2.1159,2.9105] and [1.6777,3.1455] at m = 2, 
and [1.1011,3.0452] at m = 3, whereas in Fig. [6(b), we obtain for the first intricate mode 
at V! > v 2 : [2.5865,1.8154] and [2.2191,2.3176] at m = 1, and [1.4665,2.0952] at m = 2. 
The values of /3tp are plotted versus v 2 in Fig. [H] for the TM and TE modes in a grounded 
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(b) 

FIG. 8: /?tp versus v 2 for (a) TM and (b) TE modes of a grounded DNG/DPS bilayer slab of 
relative material constants like in Fig. [2j 

DNG/DPS bilayer slab of relative material constants like in Fig. [2j One can see that /?tp has 
a stronger increase at smaller values of v 2 , whereas at greater values of v 2 , it is practically 
independent of v 2 . Thus, by coating a thin DPS layer on a grounded single-layer DNG slab, 
the TP is shifted towards greater values of and v± for each TM or TE mode of order 
rn^O, the shift being stronger for the TM modes in our numerical example. 

VI. CONCLUSION 

In this paper we analyzed the TM and TE modes in a grounded slab containing DNG /DPS 
bilayers. Simple normalized dispersion relations were given for the guided and evanescent 
modes. Relations corresponding to the grounded single-layer DNG slab were refound as 
specific cases. Numerical examples were given showing dispersion curves of the lower order 
modes and the total normalized power carried on the propagation direction by the respective 
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modes. Examples of electromagnetic fields inside the grounded DNG/DPS bilayer slabs 
were also given. Snapshots obtained by the FDTD method 23|. Implicit relations for the 
normalized parameters at the TP were given for the grounded DNG/DPS and DNG slabs. 
We showed that a given mode could have several TPs, as for example, in Fig. E] there are 
modes with three TPs, the plots of the total normalized power carried by each mode on the 
propagation direction being very useful in the designation of the TPs. Interesting results 
were obtained for the grounded DNG/DPS bilayer slabs of material constants (ei,/xi) = 
(—62,-/12) [see Fig. [5] and relation (TTT)) ] and of material constants (ei,//i) = (— //2, — £2) 
(see Figs. Eland [7]). For simplicity, we considered e\\i\ > e 2 fi2 and several examples were 
given, but relations presented in terms of normalized parameters in this paper could be 
applied also to other combinations for the relative material constants of the DPS and DNG 
materials. The analysis in this paper was restricted to the simplest case of lossless material 
constants. When absorption is necessarily taken into account, then the predicted behavior 
can qualitatively change. 
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